Systematic analysis of symmetry energy effects in the neutron star crust properties 
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The functional form of the nuclear symmetry energy in the whole range of densities relevant for 
the neutron stars is still unknown. Discrepancies concern both the low as well as the high density 
behaviour of this function. By use of Bezier curves three different families of the symmetry energy 
shapes, relevant for different density range were introduced. Their consequences for the crustal 
properties of neutron stars are presented. 
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I. INTRODUCTION 

The basic quantity in the description of infinite nuclear 
matter filling out the interior of neutron star is the energy 
per particle expressed in term of baryon number density 
n n and isospin asymmetry a = — " p of the 



it — iL p 
system: 



E(n,a) = V(n) + E s (n)a 2 + C(a 4 ) (1) 



Instead of a it is useful to use proton fraction x, and 
then a = (1 — 2x). Here we assumed the only con- 
stituents of stellar matter are nucleons and leptons: elec- 
trons and muons. Around and above the nuclear density 
uq = 0.16 fm nucleons and leptons form a quantum 
liquid, which stands for liquid core of the neutron star. 
Slightly below no matter cannot exist as a homogeneous 
fluid - the one-phase system is unstable and the coexis- 
tence of two phases is required. At these densities matter 
clustcrizes into positive nuclei immersed in a quasi-free 
gas of neutrons and electrons. Model calculations show 
this type of matter forms a Coulomb lattice with solid 
state properties and corresponds to the crust covering 
liquid core of a star. The presence of the crust is af- 
firmed by the glitching phenomenon observed for some 
pulsars. For typical NS masses, between 1-2 M Q , the 
most of the stellar matter is occupied by the core, so the 
global parameters like the mass, radius, moment of iner- 
tia are completely determined by the functional form of 
the Eq. (HJ). Whereas the isoscalar part V(n) corresponds 
mainly for the stiffness of Equation of State (EOS) which 
is relevant for the maximum mass of NS, the isovector 
part E s (n) is responsible for the chemical composition of 
the matter. Through the /3-equilibrium equations 



4(1 - 2x)E s (n) = n n - Hp = u e = 



(2) 



the proportions of all particles are determined. The sym- 
metry energy is also relevant for the crust-core transition 
in NS as it was shown in [l[. It's role can be explicitly 
seen if one looks at the compressibility under constant 
chemical potential relevant for the stability of homoge- 
neous beta equilibrated nuclear matter: 



= n 2 (£7> 2 + V") + 2n(E' s a 2 + V) - 



2a 2 E l2 n 2 



(3) 



and when > the matter is stable. This formula 
accounts for the bulk approximation and can be improved 



by inclusion of finite size effects like Coulomb and surface 
contributions. Such corrections were studied in Q were 
it was shown that 



v{Q) 



v + 2(47re 2 /3) 1/2 - j3k% F 



(4) 



is the minimal value for stable density modulations for 
the Q momentum. Stability of matter is given by the 
condition v(Q) > 0. Note that the vanishing of vq is 
equivalent to the vanishing of since they are related 
by 



v (n) 



8K u (n)E s (n) ( dfin 



(5) 



Both of these approaches consider stability of a one phase 
system against density fluctuations. However there exists 
another approach based on treating the NS crust as a 
two component system subject to the Gibbs conditions 
for mechanical and chemical stability expressed by: 



P =P 



ii 



He = M 



II 



(6) 



Where the first component (I) corresponds to clusters 
composed of protons and neutrons immersed in a pure 
neutron liquid (II) component. Both phases are per- 
meated by degenerated electrons. As density increases 
towards the star interior the two phase system can no 
longer exist and it signals the crust-core transition. These 
approaches correspond to three different critical densities 
n c (Kp), n c (Q), and n c (l <-> 2) which will be presented 
for various models in the following sections. 

In this work we are interested in the crustal proper- 
ties, so it seems natural to focus on the symmetry energy 
form around saturation density no, to which the critical 
density n c is closely located. However, as was suggested 
in Q the crust is affected by the star compactness which 
depends on the E s form at densities much higher than 
no. Furthermore recent experimental measurements Q 
show non-standard behaviour of E s at densities much be- 
low no. All these issues lead us to the idea to somehow 
"factorize" the shape of E s to see how symmetry energy 
at different ranges of density affects the NS crust. It ap- 
peared to be possible to define models in which the E s 
is changing at chosen range of density whereas the rest 
of its shape is kept the same. We have constructed the 
three main families of models with varying shape at very 
low, n — > 0, intermediate n w uq and very high densi- 
ties n > The work is organized as follows: in the 



2 



section II we present recent experimental data useful in 
constraining the form of the symmetry energy, in the sec- 
tion III the construction of different models is shown. In 
the last two sections the results of the calculation and 
the astrophysical constraints coming form NS observa- 
tions are discussed. 



II. MEASURED SYMMETRY ENERGY 
PROPERTIES IN LABORATORY 
EXPERIMENTS 

The symmetry energy E s properties around saturation 
density can be inferred from laboratory experiments. For 
nuclei, it corresponds to the volumetric symmetry energy 
term S v in the liquid droplet model [5]: E s (tlq) = S v , 
which also contains the surface contribution to the sym- 
metry energy S s . By measuring in nuclei these two quan- 
tities is possible to derive E s (no) and L which is related 
to its slope 



3n 



dE s (n) 



dn 



(7) 



Not precisely available at the moment from experiments 
but also relevant for NS matter is the curvature of the 
symmetry energy: 



9n 



dn 2 



(8) 



Neutron rich nuclei feature halos composed mainly of 
neutrons while most of the protons stay in nuclear cores. 
Such neutron skin thickness is dependent on both volu- 
metric and surface symmetry energies (by means of the 
coefficient S v /S s ) and provides another way of measuring 
E s (jiq). S v itself can be extracted from heavy ion colli- 
sions where isospin diffusion occurs: two colliding nuclei 
with opposite neutron abundances exchange components 
reaching isospin equilibration. Both, neutron skin thick- 
ness measurements [f| and isospin diffusion Q shows that 
L is placed in the range between 40 and 80 MeV. 

Experiments allows for measurements of asymmetric 
matter compressibility K asy being related to K s by the 
approximate relation K s « K asy + 6L. Recent experi- 
mental analysis [8l4icj| shows that K asy takes values be- 
tween -650 and -400 MeV. It means that, including the 
discrepancies of L, one may estimate that K s is in a broad 
range between -400 and 100 MeV. 

The density around no is the natural place where the 
efforts in determining of symmetry energy behavior are 
focused. However even if the values of slope L or the 
curvature K s achieve satisfactory accuracy the symmetry 
energy shape remains uncertain in both very low and 
very high density which are relevant for NS properties. 
For some years, the very low density part of the nuclear 
energy has been studied as well 0, EH • In these works 
the symmetry energy is probed at n « 10 -2 fin -3 and 
appears to take large values ~ 10 MeV at this range of 



density. The result is very interesting. First, it allows for 
going beyond the poor and not very restrictive expansion 
around the no, and secondly it is highly relevant for the 
crust-core transition in neutron star. 



III. BEZIER CURVES FOR THE SYMMETRY 
ENERGY 

For a good description of the neutron star properties 
it is important to know the symmetry energy functional 
form in the whole range of densities. From the aforemen- 
tioned measurements one can pin down the symmetry 
energy properties around saturation but both low and 
high density parts are not well determined. Therefore 
there is large freedom to assume the behavior of E s (n) 
and introduce different parametrization. 

Previously used parametrization (often polynomials) 
were not convenient since once the same saturation point 
properties are fulfilled the values far from saturation are 
fixed, thus sometimes leading to unphysical results, like 
unstable EOS (MDI parametrization) for large negative 
vales of E s at high densities @ . Polynomial interpolation 
presents some advantages as it may produce any shape 
but might present problems when derivatives are com- 
puted. At the boundary of the domains of interpolation 
the analicity is lost. In particular for pressure and com- 
pressibility the first and second derivatives are required 
and such polynomial interpolation leads to artifacts in 
the EOS. On the contrary Bezier curves [l2| allow for 
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FIG. 1. Bezier Curve (thick blue) for seven control points 
(brown). The thin lines represent symmetry energies for the 
MDI models shown as a reference. 

construction of different E s shapes that respect satura- 
tion properties, like E s (uq) and L, but can take any value 
at low and high density values, see Fig. [T] The resulting 
Bezier curve is then an analytical function that offers big 
advantages over interpolated functions. 

Bezier curve of degree n is based on n+ 1 control points 
Pi, where i = 0, 1, 2 . . . n. For implementation purposes, 
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an arbitrary number of control points can be chosen to 
cast the curve shape as desired by means of the following 
relation: 



b(*)=E " a-*) 



i=0 



VP,- 



t G [0, 1] (9) 



where ("J is the binomial coefficient. 

The constructed models for E s are presented in the fol- 
lowing subsections and their control points are tabelar- 
ized in the appendix. As we are interested only in the 
symmetry energy effects, the isoscalar part V(n) in the 
Eq. (fTJ) , required for the of the full nuclear model, is kept 
all the time the same. We used the functional form taken 
from [l3l ] for which the symmetric matter compressibility 
is K = 240 MeV. 



A. Low density symmetry energy effects and the 
neutron clusterization problem 

In an earlier work the authors implemented a set of 
models motivated by the low density behaviour of the 
symmetry energy that were introduced to explore NS 
properties . However a detailed analysis of phase tran- 
sitions revealed pathological properties. Fig. [5] shows the 
phase diagram of a model similar to the earlier klO in 
which E s (n ~ 0) = 10 MeV and for which the thermo- 
dynamical Gibbs conditions hold. In it, the region to the 
left of the spinodal line (in red) where the proton frac- 
tion is zero (pure neutron matter) is reached by isobaric 
lines of negative pressure which indicates clusterization 
of neutrons. This is in disagreement with common knowl- 
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FIG. 2. AE-kW model in which E s (0) = 10 MeV. Clus- 
terization of pure neutron matter occurs when the energy per 
baryon is given by E(n, x). 

edge that neutrons never clusterize. Here we propose a 
method avoiding this pathology. In order to keep finite 
values of E s at vanishing density it is then necessary to 



correct the whole expression for the total energy of nu- 
clear matter. First, we introduce the following symmetry 
energy expression: 



E s {n) 



E PALu^ 

Ef ALu {n) 



E. 



Bez 



• k (n) if n < no 
if n > n , 



where E^ ALu (n) is the symmetry energy introduced 
in the work [13j with the interaction part F{u) = u. 
Ef ezier > k (n) is a Bezier curve that modifies the low den- 
sity part in such a way that Ef Szier ' k (0) + E^ ALu (0) = k, 
k being any value in MeV. In this way, E s (n) can take 
the k = E s (0) value for which the corresponding Bezier 
curve is designed for. With it, the energy per particle 
becomes: 



E(n, x) = V(n) + (1 - 2x) 2 E s (n). 



(10) 



The above formula still needs to be corrected for neutron 
clusterization. In order to avoid neutron clusterization, 
the following correction is defined: 



AE{n) = 



j^jBezier 



(n) if n < no 
if n > n , 



so that the final form of the energy is given by 
E(n, x) = E{n, x) - AE(n). 



(11) 



We call AE-k models when the above modifications are 
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FIG. 3. Energy per particle E(n,x) for the AE-klO model 
with correction (solid lines) and without it: E(n, x) (dashed 
lines). The blue lines correspond to pure neutron matter while 
the purple ones correspond to symmetric nuclear matter. 

applied. In fact, the symmetry energy can be recovered 
by definition: 

E s {n)=E(n,x = Q)-E(n,x = -). (12) 

and captures the features of the k models (finite sym- 
metry energy values at almost zero baryon density). At 
saturation point these models have the same values as 
the PALu model (in MeV): 



E s {n ) =30, 



L = 77, 
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B. Different symmetry energy slopes at saturation 
density 

The importance of the slope of the symmetry energy 
for the NS crust-core transition point can be easily seen in 
Eq. ([3]) . The last term in that equation is the one respon- 
sible for the stability breaking and includes the square of 
the symmetry energy slope. The measured value of L 
will influence the behaviour of E s away from saturation. 
In this section we introduce two families of models with 
the same symmetry energy forms at high and low den- 
sities but with variable slope at saturation point. There 
at saturation, they share the common values 

E s (n ) = 31, K s = 0. 



FIG. 4. AE-klO model in which E B (0) = 10 MeV. Clusteriza- 
tion of pure neutron matter is now impossible for this model 
when the AE correction is included. Close to the spinodal 
(in red), the system splits into two phases, (1) - pure neutron 
(2) - nuclear matter represented by dotted green lines. 



Let us recall that the Ef ezter curves where defined in 
such a way to preserve the known properties of symmet- 
ric nuclear matter, i.e., its compressibility at saturation 
and its minimum. These conditions are fulfilled because 
the following derivatives with respect to baryon number 
density vanish at saturation: 



whereas the slope L takes values in the range between 
40-120 MeV corresponding to those reported by various 
experiments [2l| . These properties are established by 
joining two Bezier curves describing the low and high 
density parts joint at saturation point with continuous 
derivatives up to second order. At low densities all of 
them go to zero while at high densities they follow the 
behavior of either b or c MDI models Q ■ The two families 
are presented in Fig. |5J The reason of introducing the 
b or c-like behaviour at high densities comes from the 
analysis of direct URCA constraint which is discussed in 
detail in the subsection D. 



AE'(n ) = 0, 
AE"(n Q ) = 0. 



(13) 
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FIG. 5. Nuclear symmetry energy for AE-k models. The 
curves correspond to E s (0) = 0, 4, 8, 12 MeV. 



FIG. 6. Different symmetry energy shapes for the L-high-c 
models (top) and for L-high-b (bottom). 
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C. High density region symmetry energy effects 

The critical density for the crust-core transition is the 
most relevant quantity for determination of the NS crust 
properties. Its value determines the crust thickness and 
moment of inertia. The crust-core position depends only 
on the behavior of the symmetry energy below satura- 
tion. However the crust thickness is not only determined 
by the critical density but also by the global parameters 
of the NS like its mass and radius. Those global param- 
eters depend on the high density part of E s . To explore 
such effects we construct a set of models with the same 
low density part implying the same crust-core transition 
but with different high density part. In relation to that, 
the authors of [l4| derive an approximate formula assum- 
ing a polytropic equation of state for a thin, light crust 
but deviations from those assumptions can lead to com- 
pletely different results and are due to the high density 
E s region, as will be presented in the following sections. 
We introduce the high-E s models which have the follow- 
ing values (in MeV): 

E s (n ) = 31, L = 60, K s = 

and are composed of two Bezier curves describing both 
the low and high density parts joint at saturation point 
with continuous derivatives, just like in the previous 
models for variable slopes. They are named b, c, d since 
they possess the same high density values as the corre- 
sponding b, c, d MDI models [2] . The e model is an ex- 
treme case whose values grow much faster with increasing 
density than the rest. Fig. [7] shows E s for all of them. 
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FIG. 7. Symmetry energies for the high-E s models. 



IV. ASTROPHYSICAL OBSERVATIONS AND 
SYMMETRY ENERGY CONSTRAINTS 

Neutron star models must be compared with astro- 
physical observations to prove their validity. Although 
there are many astrophysical processes, not all them al- 
low for precise information about the NS properties, like 
simultaneous accurate mass and radius measurements. 
However, recently the highest neutron star mass has been 



precisely measured [l5| with a value around 2 M Q plac- 
ing a strong constraint to many of the existent equations 
of state. Even though E s contributes to the determi- 
nation of the maximum NS mass for an EOS it is the 
isoscalar part V that plays the major role. Therefore this 
measurement is not so stringent for E s . Another perti- 
nent result regarding NS crusts is the upper bound on 
the crustal moment of inertia for the Vela pulsar derived 
from an analysis performed by Link in telling that 
such quantity should be higher than the 1.4% of total 
moment of inertia of this star. Different E a forms will 
result in different crust properties to be tested against 
this condition. As a final constraint for the models pre- 
sented here the Direct Urea cooling is considered. Low 
mass NS do not cool by the DUrca process according 
to [20j . The proton fraction of the star x should not go 
above the DUrca proton fraction threshold xdu for those 
low NS masses. The DUrca threshold xdu is weakly de- 
pendent density function and takes values between 0.11 
and 0.15. From the /^-equilibrium equations, Eq. fol- 
low values of x and neglecting muons one may get the 
following estimation for the proton fraction in matter 



6 + n/n (83 MeV/E s ) 3 ' v ' 

This expression implies that one natural way to fulfill this 
constraint is restricting E s to low values, smaller than 80 
MeV, so the resulting x always stay below xdu- That 
is the reason we introduce the L-high-b models - the sec- 
ond sub-class in the L-high family, which presents soft 
behaviour at high density. As we will see in the following 
sections, this possibility has important implications for 
the EOS when satisfying the maximum NS mass mea- 
surement. 



V. NEUTRON STAR PROPERTIES 

In order to construct the neutron star profile we need to 
include an EOS in the whole range of densities from zero 
to high densities characteristic of the very central part of 
the star. In this study we have considered a liquid core 
described by EOS based on the symmetry energy forms 
presented in the previous section with a common isoscalar 
part given by the PAL parametrization. As for the crust, 
we have taken the SLy EOS composed of different parts 
whose table can be found in [l6j]. These two EOSs are 
joint at the point of equal pressure and density. This 
construction is not thermodynamically consistent since 
there is a discontinuity in chemical potentials but can 
be treated as an approximation that does not influence 
the overall macroscopic NS properties. In our approach 
the bottom edge of the crust is determined by the crit- 
ical density n c {K^) which is very close to n c (l O 2). 
The NS total moment of inertia and that of its crust are 
calculated by the following expressions, derived in the 
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framework of General Relativity |l9[, as: 



I = 



J 



1 + 2GJ/R 3 c 2 ' 



1 - 2GI/R 3 c 2 
1 - 2GM/Rc 2 



^■Icrust — -z(M crus tR 2 ) 



(15) 



the total mass and the mass of the core, as determined 
by n c . Interestingly, the consequences of choosing n c (Q) 
are systematic effects of lowering the NS crust thickness 
and its moment of inertia. The effect was shown in 1171. 



A. AE-k models 

In these models the E s forms take finite values for very 
low densities in the range between — 12 MeV sharing 
the same high density behavior and same L value. They 
result into an acceptable maximum mass (of about 2 Mq) 
and radii between 12 — 13 kilometers. Table Q] shows the 
transition densities for each k value: as k grows the tran- 
sition density n c lowers. A study of the behavior of the 
compressibility curves in this set of models show that 
as we increase the values of the symmetry energy at low 
densities, the point of vanishing occurs at lower and 
lower densities to finally does not appear for k greater 
than 14 MeV. This results into thinner and thinner neu- 
tron star crust and could be extended to crustless neutron 
stars, which are, however, not expected to exist accord- 
ing to observations. The crustal properties appear to 



TABLE I. Crust-core transition densities (fm 3 ) for the AE-k 
models. 

model E s (0) n c (Q) n e (_K" M ) n e (l o- 2) 



AE-kOO 
AE-W2 
AE-WA 
AE-k06 
AE-k08 
AE-klO 
AE-kl2 
AE-kU 




2 
4 
6 
8 
10 
12 
14 



0.0816675 
0.0690383 
0.053707 
0.0399434 
0.0306162 
0.0248644 
0.0209877 
0.0173067 



0.0930797 
0.0814262 
0.0647734 
0.0473238 
0.0351069 
0.0278587 
0.0233639 
0.0198872 



0.0942588 
0.0820099 
0.0647989 
0.0474521 
0.0353349 
0.0288262 
0.0257871 
0.0215231 



be very sensitive to the varying portion of the symmetry 
energy shape, specially the crustal moment of inertia, 
which is shown in Fig. [8] The most likely value for the 
low density E s part derived from experiment is around 10 
MeV. If it is true, it will point to a very low mass of the 
Vela pulsar, much below 1 M Q , according to the models 
here. Such low masses are not favored by the supernova 
explosion scenarios which produce new born neutron star 
with masses around 1.5 Mq. This result is similar to the 
already reported in Q but there it could be interpreted 
as an inconsistency of the model coming from the clus- 
terization of pure neutron matter. For the models here 
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FIG. 8. Results for the AE-k models. Top: crust thickness. 
Bottom: fractional moment of inertia in the crust where the 
red line represent the lower bound for crustal / of the Vela 
pulsar. Each curve corresponds to a specific value of k — 
E s (0). 



this inconsistency has been removed but the same effect 
persists. Summarizing, the above results signal that the 
large values of E s at n — > are in contradiction with 
neutron star observations. 



B. L models 

A first look at the E s form makes an impression that 
differences between different models in the L-high-c fam- 
ily are not large since they almost overlap and coincide 
at low and high density regions. Basically the only dif- 
ference between them is the slope at saturation point. 
From table [TT| we see that the critical density for crust- 
core transitions is not clearly correlated to L values: for 
both high and low L the critical density is large whereas 
for the intermediate L is lower. The behavior of the 



critical density with L was already analyzed in [22j and 
shows differences when going beyond the parabolic ap- 
proximation in the nuclear energy used here. In our ap- 
proach critical density differences are not very large but 
on the contrary the NS parameters are highly sensitive 
to the values of L. Both global and crustal properties of 
NSs are affected by such the different values of L. One 
may see that the star compactness j3 = GM/(Rc 2 ) is es- 
sentially changed by L specially for low NS masses. It 
could be naturally expected since in low mass NS central 
densities are relatively low and the low density E s values 



7 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

Mass[solar mass] Mass[solar mass] 



FIG. 9. Neutron star features for the L-high-c models. 



TABLE II. Crust-core transition densities for the L-high-c and 
L-high-b models. 

model n c (Q) n e (i"C M ) n c (l -H> 2) 

L40-high-c 0.101234 0.107808 0.113369 

L60-high-c 0.0918645 0.100315 0.103118 

L90-high-c 0.0876315 0.100575 0.101674 

L120-high-c 0.111214 0.139422 0.143345 

L40-high-b 0.101234 0.107808 0.113369 

L60-high-b 0.0918645 0.100315 0.103118 

L90-high-b 0.0876315 0.100575 0.101674 

L120-high-b 0.111214 0.139422 0.143345 



mainly determined by L have to be relevant. The E s 
slope contributes to the pressure so increasing L makes 
the EOS stiffer and decreases the compactness. On the 
contrary, a soft equation of state amplifies the effect of 
gravity on the star as a result of a more compact star. It 
is interesting that even if the crust-core transition does 
not behave monotonically with L, the crustal properties 
like thickness and moment of inertia depend clearly on 
the E s slope. The larger L is, the thicker the NS crust 
becomes. The effect is more pronounced for the moment 
of inertia carried by the crust. Small values of L corre- 
spond to softer EOS resulting in a more compact star. A 
more compact object has stronger gravity in the crustal 
region, therefore the net effect is that the crust is more 
squeezed and contributes less to the total moment of iner- 
tia. Summarizing, in the case of L models, the crust-core 



transition point n c does not play a major role in determi- 
nation of the NS crustal properties but rather the gravity 
controlled by L is the main factor here. 



C. high-E s models 

In this section we present results for the third family. 
All models in this family share the same shape of the 
symmetry energy up to saturation point, which means 
that they show the same critical density for the crust- 
core transition point, shown in Table [TTTT Therefore the 
crustal properties are only affected by the high density 
shape of E s . As one may see from Fig. [TUl the compact- 
ness changes a lot for massive stars in the contrast with 
the previous family where this property changed for low 
massive stars mainly. The form of the symmetry energy 
at high densities changes the compactness of the NS so 
that the influence of gravity is different in the crust re- 
gion in each model. For larger compactness the crust is 
more squeezed, has lower thickness and contributes less 
to the total moment of inertia, what is shown in Fig. [TT] 
Those differences are more or less constant and do not 
depend on the total mass of a given NS. For moments 



TABLE III. Crust-core transition densities for the high-E s 
models. 

model n c (Q) n c (i(' M ) n c (l -H> 2) 

high-E s models 0.0918645 0.100315 0.103118 
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FIG. 10. Mass vs radius relation for the high-E s models. 

of inertia the differences are larger then for the crust 
thickness. Here we must emphasize that those differences 
comes only from the portion of E s at very high density 
which is responsible for the different compactness of the 
star. The influence of the compactness on crustal prop- 
erties of NS was estimated by an approximate formula 
for the moment of inertia derived in [14| : 

AI crust ^ 28nP t R 3 (l- 1.67/3 - 0.6/3 2 ) 
I ~ 3Mc 2 

*(i+ 2Ft(i+w 

\ n t m b e L P z J 

where P t = p(n t ) and n t (called n c in this work) are 
the values of pressure and baryon number density at the 
crust-core transition and P is the compactness param- 
eter. In the Fig. [TT] we compare the results of this ap- 
proximation and we see discrepancies that show its range 
of applicability. It means that the scaling of the crustal 
properties with the compactness is not simple and this 
formula must be taken with care specially in case of thick 
crusts. 



D. Direct Urea constraint 

According to [20| neutron stars with masses below 1.35 
M Q should not cool by the direct Urea process. The pro- 
ton fraction of the star x should not go above the DUrca 
proton fraction threshold xdu for those low NS masses. 
It is mainly determined by the symmetry energy, imply- 
ing that one way to avoid violating this constraint is re- 
stricting E s to low values so the resulting x always stay 
below xdu (see [23[ for detailed discussion). That in- 
deed does not happen for most of the models studied in 
this work. The energy expression, Eq. (fT]), whose sym- 
metric part is based on the PAL parametrization leads 
to rather soft EOS, so the central density may be eas- 
ily attained by low massive stars. The critical masses 
satisfying the DUrca constraint are: 1.2 M Q for AE-k 
models, 0.5 - 1.0 M Q for L-high-c models and 0.5 - 1.3 
M Q for high-E s models. In the particular case of the 
L models the higher density part of E s is meant to be 
completely arbitrary and surely one can incorporate the 
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FIG. 11. Effects of the different symmetry energy forms in the 
high-E 3 models for crustal thickness (top) and moment of in- 
ertia (bottom). The dashed lines correspond to the AI cru3t /I 
derived by the approximate formula Eq. (|16p 




0.0 0.1 0.2 0.3 0.4 



n[fm" 3 ] 

FIG. 12. The proton fraction and DUrca constraint for the 
L-high-b family. 



cooling limitation easily and improve them, by means of 
Bezier curves. Figure [12] shows the proton fraction for 
the L-high-b models whose high density part of the sym- 
metry energy behaves like the b MDI model. There it 
can be seen that the DUrca cooling never sets in since 
the proton fraction always stays below the DUrca thresh- 
old. This constraint therefore favors E s forms which do 
not grow quickly with increasing density and stay bound 
from above. However such soft E s leads to the maximum 
mass not greater than 1.85 Mq which is in conflict with 
observations. 

Nevertheless, there exists another possibility apart 
from modifications to E s to fulfill the aforementioned re- 
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striction. The symmetry energy contributes to the stiff- 
ness of an EOS, but the main contribution comes from 
the isoscalar part of the nuclear energy. As already men- 
tioned, in most of the models with arbitrary E s used here 
the DUrca cooling sets in for unacceptable masses, due 
to the particular form of E s . But if it turns out that E s 
is such that x never goes above the DUrca threshold x du 
(like in the L-high-b models) then it's most likely that the 
isoscalar part of the PAL parametrization be wrong. It 
should be then properly replaced by an expression that 
enhances the stiffness of the EOS that would not only sat- 
isfy the DUrca constraint but at the same time contribute 
to the appearance of high mass neutron stars. It is worth 
mentioning that even if the PAL formula is highly biased, 
it has been useful to spot the E s contribution to the NS 
physics. On the contrary, if a PAL-like isoscalar part 
is correct, then the symmetry energy can be very much 
constrained by both the cooling and the maximum mass, 
resulting in moderate values at low densities (to satisfy 
DUrca) and then increasing abruptly so high mass NS 
are created within the model. This result could be in 
fact very stringent, and its validity is to be studied and 
confirmed with the upcoming observations. As a final re- 
mark it is important to mention that future approaches 
to the EOS should take the cooling phenomenon as an 
important test. 

VI. CONCLUSIONS 

In this work we were interested in the role of the sym- 
metry energy form played in the NS crust properties. 
The crust-core transition occurs well below the saturation 
density, so the present knowledge of the E s parameters 
at this density is not sufficient to determine the transi- 
tion point exactly. Moreover the symmetry energy form 
at higher densities, although not affecting the crust-core 
transition point directly is also relevant. One may say the 
whole shape of the E s has to be taken into account. In 
order to extract the influence of different portions of the 
E s shape we have constructed the three different families 
of symmetry energy parametrization corresponding the 
three range of densities: well below the saturation point, 
around n and highly above. In each family one chosen 
feature was changing whereas the remaining portion of 
E s (n) was kept the same. This was achieved by use of 
Bezier curves parametrization. 

In the low density regime (A-E'-fc-models) the charac- 
teristic feature was the value of E s (n — > 0) being not as 
it is suggested by some recent measurement. It appeared 
that high values of E s at low density are questionable as 
they lead to a very thin crust, difficult to be reconciled 
with observations of the Vela pulsar. 

In the intermediate range of density L-models the key 
quantity was the slope of symmetry energy L taken in 
the range from 40 to 120 MeV. The crust thickness and 
moment of inertia is very sensitive and highly increases 
with the value of L. 

For the high densities we constructed the models pre- 



serving the same crust-core transition density but pre- 
senting completely different behaviour well above uq, 
they were called high-E s models. It appeared that values 
of E s at high density affect the crustal properties essen- 
tially. Such effect comes from different compactness of 
a star. However it is difficult to find a simple scaling of 
crust properties with M / R, especially when the crust is 
thick. 

The direct Urea constraint for low mass NS is related 
to the symmetry energy since it serves to determine the 
proton fraction inside the star. If the high density behav- 
ior of E s is bounded from above to satisfy this constraint 
a 2 M NS cannot be created by the models used here. A 
stiffcr isoscalar part V of the nuclear energy per particle 
is then necessary to produce higher masses. Therefore 
both the heaviest observed NS and the DUrca cooling 
condition allow for constraining the EOS, since they are 
to be satisfied simultaneously. Solving this issue is surely 
an interesting implementation for a future work. 

During the preparation of the text the authors have 
found the work by Lattimer et al. (24j which combines 
the majority of laboratory measurements. It is concluded 
there that the most reliable values of S v « 32 MeV, L ks 
50 MeV with an error of a few MeV. This means that in- 
dependence on the crustal properties is almost removed. 
In our work, to avoid additional complexity, we have fixed 
the second derivative K s = for most presented models. 
Now, when the first derivative is pinned down it seems 
interesting to explore different K s values, what is planned 
for the future. 
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Appendix: Bezier points for various models 



TABLE IV. Bezier control points for the AE — k models. 
AE - kOO is simply the PALu model. 

model Pp Pi P 2 P 3 P 4 

AE-k02 (0, 2) (0.04, 1) (0.08, 0) (0.12, 0) (0.16, 0) 

AE -k04 (0, 4) (0.04, 2) (0.08, 0) (0.12, 0) (0.16, 0) 

AE - k06 (0, 6) (0.04, 3) (0.08, 0) (0.12, 0) (0.16, 0) 

AE - k08 (0, 8) (0.04, 4) (0.08, 0) (0.12, 0) (0.16, 0) 

AE-kW (0, 10) (0.04, 5) (0.08, 0) (0.12, 0) (0.16, 0) 

AE-kl2 (0, 12) (0.04, 6) (0.08, 0) (0.12, 0) (0.16, 0) 

AE - kU (0, 14) (0.04, 7) (0.08, 0) (0.12, 0) (0.16, 0) 
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TABLE V. Bezier control points for the L models which are composed of two Bezier curves joint at no. 



model 


Po 


Pi 


P 2 


P3 P4 Ps 


L40-high-c 


(0, 0) 


(0.0528, 22.066) (0.106, 26.466) 


(0.16, 31) 


(low density) 










L40-high-c 


(0.16, 31) 


(0.24, 37.666) 


(0.32, 44.333) 


(0.8, 235.372) (1.28, 455.583) (1.92, 819.168) 


(high density) 










LoU-nign-c 


(0, 0) 


(0.0528, 17.6) 


(0.106, 24.2) 


(0.16, 31 J 




(0 16 31 1 


(0.24, 41) 


(0.32, 51) 


(0 8 235 372"! fl 28 455 5841 Cl 92 819 1681 


t on Viin-Vi r> 
ijyu-nign-c 


In n\ 

(U, U) 


(0.0523, 10.9) 


(0.106, 20.8) 


l^U. ID, ol ) 




(0 16 31 1 


(0.24, 469) 


(0.32, 61) 


f0 8 235 3721 fl 28 455 5841 fl Q2 81 Q 1 681 


ijizu-nign-c 


In n\ 

(U,U) 


(0.053, 4.2) 


(0.106, 17.4) 






(■0 16 31 "1 


(0.24, 51) 


(0.32, 71) 


(OH 235 3721 fl 28 455 5841 fl Q2 81 Q 1681 


L40-high-b 


(0, 0) 


(0.0528, 22.066) (0.106, 26.466) 


(0.16, 31) 




(0.16, 31) 


(0.24, 37.666) 


(0.32, 44.333) 


(0.8, 72.239) (1.28, 79.064) (1.92, 71.955) 


T c\c\ Viio-Vi h 

JL-UU-Illg,-! - U 


("n nl 


(0.0528, 17.6) 


(0.106, 24.2) 


fn ifi qi"! 




(0.16, 31) 


(0.24, 41.000) 


(0.32, 51.000) 


(0.8, 72.239) (1.28, 79.064) (1.92, 71.955) 


L90-high-b 


(0, 0) 


(0.0523, 10.9) 


(0.106, 20.8) 


(0.16, 31) 




(0.16, 31) 


(0.24, 46.000) 


(0.32, 61) 


(0.8,72.239 ) (1.28, 79.064) (1.92, 71.955) 


L120-high-b 


(0,0) 


(0.053, 4.2) 


(0.106, 17.4) 


(0.16, 31) 




(0.16, 31) 


(0.24, 51) 


(0.32, 71) 


(0.8, 72.239) (1.28, 79.064) (1.92, 71.955) 



TABLE VI. Bezier control points for the high — E a models. 



model 


Po 


Pi 


P 2 


P 3 P 4 Ps 


low density 


(0, 0) 


(0.0528, 17.6) (0.106, 24.2) 


(0.16, 31) 


high density 


L60-high-b 
L60-high-c 
L60-high-d 
L60-high-e 


(0.16, 31) 
(0.16, 31) 
(0.16, 31) 
(0.16, 31) 


(0.24, 41) 
(0.24, 41) 
(0.24, 41) 
(0.24, 41) 


(0.32, 51) 
(0.32, 51) 
(0.32, 51) 
(0.32, 51) 


(0.8, 89.12) (1.28, 82.86) (1.6, 71.92) 
(0.8, 207.26) (1.28, 440.014) (1.6, 626.88) 
(0.8, 301.097) (1.28, 813.098) (1.6, 1146.13) 
(0.8, 547.93) (1.28, 1244.98) (1.6, 1673.62) 
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